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^ \ Abstract 

^ . It has been pointed out that the value of the gravitational "constant" in the early 

QQ ' universe may be different from that at present. In that case, it was conjectured that 

\ primordial black holes may "remember" the value of the gravitational constant in the 

■ early universe. The present analysis shows that this is not the case, at least in certain 

' contexts. 

o 
o 



1 Introduction 

In the early universe, black holes may have been formed due to inhomogeneous initial conditions, phase 
transitions, or other mechanisms. They are called primordial black holes (PBHs). The mass of the PBHs 
O ' is of order the mass contained within the Hubble horizon at the formation epoch. Hawking |l[] discovered 

^ , that black holes emit radiation due to the quantum effects of the curved spacetime. PBHs therefore lose 

^ ' mass and those lighter than ~ lO^^g may contribute to the cosmic gamma-ray background and spoil the 

success of the big bang nucleosysthesis scenario. PBHs heavier than that may dominate the density of 
the universe at present. Therefore the fraction /3 of the universe which may have gone into PBHs can be 
constrained from various cosmological observations (Carr ||] ) . 
\^ ' It has been pointed out that gravity in the early universe may have been dilatonic, since this arises 

5^ , naturally as a low-energy limit of string theory. In such theories, the gravitational "constant" is given 

by a function of a scalar field which couples non-minimally with gravity. In particular, the gravitational 
constant in the early universe may have been different from that at the present epoch. The Brans-Dicke 
theory is the simplest such theory of gravity. Barrow [|j considered PBHs in Brans-Dicke theory and 
proposed two scenarios. In scenario A, the gravitational constant at the black hole event horizon is 
always the same as that at the cosmological particle horizon, i.e., GEH{t) = Gpnit). In scenario B, the 
gravitational constant at the black hole event horizon is constant with time and therefore always the 
same as that at the formation epoch, i.e., GEH{t) — GEH{tf) = GpH{tf). Scenario B corresponds to 
what is called the gravitational memory, while scenario (A) entails no gravitational memory. 

Barrow and Carr Q obtained observational constraints on (3 in both scenarios. They assumed that the 
energy emission due to the Hawking radiation is determined by the temperature Th = {?>ttG EH{t)M)^^ , 
so that the mass loss rate is M ~ —{GEH{t)M)~'^. The results for the two scenarios can be significantly 
different. Subsequently, Carr and Goymer Q argued that the true situation may be intermediate between 
Scenario A and Scenario B and this allows other possibilities. Therefore, in order to deduce the history 
of the early universe from the present observations, it is very important to examine which scenario is 
realized. 

Recently, Jacobson ^ has argued that there is no gravitational memory for a Schwarzschild black 
hole with a time-varying boundary condition which mimics the cosmological evolution of a scalar field. 
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However, this argument applies only if the black hole is much smaller than the particle horizon, so we 
must consider more general situations. The purpose of this article is to solve the problem self-consistently 
for a Friedmann-Robertson- Walker (FRW) background. 



2 Basic Equations 

The field equations in Brans-Dicke theory are given by 

Rf^i^ - ^QtJ^vR = ^^M"^ + ^ {^,tJ.<l^,y - ]^9iJ^v(t>.,a(ji'°^ + ^ (V^V,.(/) - .g^,.n0) , (1) 
V,,T^'' = 0, (2) 

where the constant lo is the Brans-Dicke parameter and (j) is the Brans-Dicke scalar field, which is related 
to the gravitational constant hy G ^ (f)^^ . In the limit uj ^ oo, the theory goes to general relativity 
with a minimally coupled scalar field. The constraint w > 3300 is required by recent observations of light 
deflection 0. However, if we consider more general scalar-tensor theories, uj also becomes a function of 
(j). In this case, lj may have been small in the early universe, even though it is large today. 

Here we calculate the time evolution of the gravitational constant by using the following post-general- 
relativistic expansion. First we assume that the scalar field is constant, i.e., (j) = 4'o- Then we find the 
general relativistic solution {gcR^TaR) which satisfies Eqs. (|l|) and (||) with G = 0q^. Next we put 
{g,T) = [gGR^Tan) in Eq. (^) and assume that dcj) — 4> — (pQ is of order uj~^ for uj 1. We can then 
obtain (j) to this order, denoting it as (j)'^^'^ i by solving the following equation: 

(4) 

where Oqji denotes the d'Alembertian in the geometry given by gcR- Substituting the solution c/)'^-' of 
Eq. into the right-hand side of Eq. (|l|), we can then determine g^^^ and T^^^ by solving Eqs. (|l|) and 
(^. Then, (^^^\ which is the solution up to 0(a;~^), is determined by solving Eq. (^) for the background 
{g^'^\ r^^)). Repeating this process, we can construct the Brans-Dicke solution from the general relativity 
solution order by order. 

Here we truncate the expansion for the scalar field at 0{ui~^). This implies that we only have to solve 
Eq. (g) and we can neglect the effect of the scalar field on the background curvature. This approximation 
was shown to be very good for the generation and propagation of gravitational waves from the collapse 
to a black hole in asymptotically flat spacetime for Brans-Dicke theory with cj > 4 H. 



3 Lemaitre-Tolman-Bondi solution 

We adopt the Lemaitre-Tolman-Bondi model as the background solution {gGRyTcR) in which the scalar 
field evolves. This model is an exact solution of general relativity which describes a spherically symmetric 
inhomogeneous universe with dust. The line element in the synchronous comoving coordinates is given 

by 

ds^ = -dt^ + A^{t,r)dr^ + R^{t,r)dn^ , (5) 

where 

dn^ ^ d0^ + sin^ Odcj)^ . (6) 

The function R(t, r) is given by 

R=[^y\t-ts{r)r/' for / = 0, (7) 
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F(r) 

^^^'> F(,) for />0, (8) 

F(r) 

^^^"^^ for /<0, (9) 

where the prime and dot denote the derivatives with respect to r and t, respectively. The function A{t, r) 
is given by 

The functions ts{r), f(r) > —1 and F(r) > are arbitrary and one of them corresponds to the gauge 
freedom. The first two relate to the big bang time and the total energy per unit mass for the shell at 
radius r. The third relates to the mass within radius r, the density of the dust being given by 

F' 

(11) 



4 Method 

We are interested in the behavior of the scalar field long after the black hole has formed. In this case, 
the characteristic method is suitable. This method was first applied to the Lemaitre-Tolman-Bondi 
background by Iguchi, Nakao and Harada p. If we introduce the retarded time u, the line element 
becomes 

ds^ = -a{u,f)du^ -2a{u,a)A{u,f)dudf + R^{u,f)dVi^ , (12) 

where a must satisfy 

d 

—a = Aa. (13) 
or 

In terms of the derivatives, d/du and d/df, the d'Alembertian is given by 

2 dip A' 1 



°^ aARdu A^r"^ ' AR 



(14) 



^ = ^W)- (15) 

We can then integrate Eq. (U) along the characteristic curves, i.e., null geodesies. With this choice of 
coordinate system, we can calculate the behavior of the scalar field long after black hole formation without 
having to impose boundary conditions on the event horizon. 



5 Models 

As we have seen, the Lemaitre-Tolman-Bondi solution has three arbitrary functions. We want to obtain 
models which describe a PBH in a flat FRW universe by fixing these functions appropriately. Here 
we adopt the following assumptions: the big bang occurs at the same time everywhere; the model 
is asymptotically flat FRW, which requires the overdense region to be compensated by a surrounding 
underdense region; the model is free from naked singularities. In order to satisfy these assumptions, we 
first choose 

ts{r) = 0, (16) 



3 



and 



— — I for r <r„ 



— exp 

r, 



(17) 



for r > r,, 



where Tc and are the curvature scale and the size of the overdense region. Using the arbitrary function 
F{r), we then set the radial coordinate r to be 



r = Rito,r) 

where t = to is the initial time. 

Before integrating Eq. (||), we have to fix the initial data for 
geneous solution given by 



(18) 



so we choose the cosmological homo- 
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— -In- 

2uj 3 <o 



(19) 



We set the initial null hypersurface as the null cone whose vertex is at {t,r) — {to,0) and regard the 
cosmological solution as the initial data on this hypersurface. Although the value of the scalar field at 
the cosmological particle horizon must be given by this solution, the value in the perturbed region may 
be different from this. We have therefore examined the effects of other initial data which are different 
from the cosmological solution in the perturbed region. However, the results are not changed much. 



6 Results 

We adopt units in which the asymptotic value of the Hubble parameter at i = is unity. We set the 
parameters which specify the initial data as r^, = 2 and r„ = 1, while the Brans-Dicke parameter is 
chosen as w = 5. It is noted that Tc must satisfy Vc ^ r^, else the overdense region is isolated from the 
rest of the universe, as Carr and Hawking jl^] pointed out. 

The results are shown in Figs. ^]^. In Figs. |l| (a) and (b), a set of initial data for the background 
geometry is shown. The initial density perturbation 6 is defined as 

Sit,r)^P-M_P^, (20) 
p(t, oo) 

In Fig. H, the trajectories of outgoing light rays in this background geometry are shown. It is shown that 
the event horizon is formed at i? ~ 10. 

The results of integrating Eq. (^) are shown in Figs, ^(a) and (b), where A(j) is defined as 

A</>^fei^. (21) 

90 

It is found that the value of the scalar field around the black hole is slightly smaller than the asymptotic 
value because of the underdensity of the surrounding region. Nevertheless, the scalar field as a whole is 
almost spatially homogeneous at each moment, in spite of the existence of the black hole. The evolution 
of the scalar field near the event horizon is well described by the cosmological homogeneous solution. 



7 Summary 

We have calculated the evolution of the Brans-Dicke scalar field in the presence of a primordial black 
hole formed in a flat FRW universe. We have found that the value of the scalar field at the event 
horizon almost maintains the cosmological value at each moment. This suggests that primordial black 
holes "forget" the value of the gravitational constant at their formation epoch. However, it should be 
stressed that this result has only been demonstrated for a dust universe in which the scalar field does not 
appreciably affect the background curvature. It remains to be seen whether the same conclusion applies 
when these assumptions are dropped. 

We are grateful to T. Nakamura for helpful discussions and useful comments. 
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Figure 1: (a) Specific energy /(r) and (b) initial density perturbation 5 at t = to- 
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Figure 2: Trajectories of outgoing light rays in the background geometry. 
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Figure 3: A0 
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<t'o)/4'o (a) along the world lines R = const and on the hypersurfaces t = const. 
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